Hitting probability for anomalous diffusion processes 
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We present the universal features of the hitting probability Q(x, L), the probability that a generic 
stochastic process starting at x and evolving in a box [0, L] hits the upper boundary L before 
hitting the lower boundary at 0. For a generic self-affine process (describing, for instance, the 
polymer translocation through a nanopore) we show that Q(x, L) — Q(x/L) and the scaling function 
Q(z) ~ as 2 — > with <j> = 9/H where H and are respectively the Hurst exponent and the 
persistence exponent of the process. This result is verified in several exact calculations including 
when the process represents the position of a particle diffusing in a disordered potential. We also 
provide numerical supports for our analytical results. 



The transfer of DNA, RNA and proteins through cell 
membranes is key to understanding several biological pro- 
cesses pp. The transport of polymer molecules across 
nanopores is also relevant in many chemical and indus- 
trial applications [3]. A fundamental question concerns 
whether a polymer, once penetrated into the pore, will 
eventually complete its transit. The answer is naturally 
formulated in terms of the translocation coordinate X(t), 
namely the length of the translocated portion of the poly- 
mer at time t [3 HI El E]. In absence of driving forces, 
the polymer dynamics is governed by thermal fluctua- 
tions. In this case, the traslocation coordinate can be 
expressed as a stochastic process X(t) that evolves in a 
box of size L (L being the polymer length) , starting from 
some initial value X(0) — x, < x < L, and terminated 
upon touching either boundary for the first time (Fig. [I] 
left). It has been shown that excluded volume effects 
hinder the polymer dynamics, so that the process X(t) 
actually undergoes subdiffusion [31 [7] . We define the hit- 
ting probability Q(x, L) as the probability of exiting the 
domain through the boundary at L, which corresponds 
to the polymer completing the translocation. 

More generally, the hitting probability Q(x,L) of a 
particle undergoing anomalous (i.e., non-Brownian) dif- 
fusion is key to understanding a variety of phenomena, 
such as the classical gambler's ruin problem in finance 
and risk management [HIS], the transport of charge car- 
riers in conductors with disordered impurities |10j and 
the breakthrough of chemical species in heterogeneous 
porous media for contaminated sites remediation [11], 
only to name a few. For ordinary Brownian diffusion, 
the hitting probability Q(x,L) = x/L is easy to com- 
pute [HMH]. The goal of this Letter is to study Q(x, L) for 
generic self-affine stochastic processes, thus going beyond 
the Brownian world. 

Since the only length scale in the problem is L, evi- 
dently Q(x,L) is a function of only the scaled variable 
xj L:Q(x, L) — Q(x/L = z). For a Brownian motion, 
Q(z) — z, a simple linear function. For a generic stochas- 
tic process X(t), Q(z) is non trivial (see for example 
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FIG. 1: Left. The evolution of a stochastic process initiated 
at X(0) — x and terminated upon exiting from the box of size 
L. Right. The function Q(z) as given by Eq. Q for different 
values of the exponent (j>. 

Fig. [I] right). The central aim of this Letter is to de- 
termine the universal features associated with Q(z) in 
the two following cases: symmetric self-affine processes, 
characterized by a power-law scaling X(t) ~ t H , with 
Hurst exponent H > 0; and a single particle diffusing in 
a disordered potential V(X). The translocation process 
belongs to the former, whereas transport in quenched 
disorder to the latter. 

It is useful to summarize our main results which are 
threefold: i) For self-affine processes, we show that gener- 
ically Q(z) ~ z^ for small z, where (j> = 9/H, and 9 is 
the so-called persistence exponent [12] of the same pro- 
cess in a semi-infinite geometry, it) For a particle diffus- 
ing in a disordered potential V(X), we provide an exact 
formula for Q(x, L) valid for arbitrary V(X) which in- 
cidentally also allows us to compute the persistence ex- 
ponent of a particle diffusing in a self-affine disordered 
potential. Hi) The function Q{z) is explicitly known for 
some anomalous diffusion processes. Amazingly, we find 
that these apparently different-looking formulae can be 
cast in the same super-universal form, when expressed in 
terms of the exponent (p. This naturally raises the ques- 
tion: how generic is this super-universality? We provide 
numerical evidences that indeed in some cases the super- 
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FIG. 2: Left. A stochastic process starting at a; leaves the 
positive half axis for the first time at tf, m denotes its max- 
imum till tf. Right. A self-affine disordered potential with 
maximum at x m : when L is large, the diffusing particle, start- 
ing at < x < L, exits the box through for x < x m and 
through L for x > x m . 



universality is violated and we discuss its limit of validity. 

Self-affine processes. To compute Q(x, L) in a box ge- 
ometry, it is useful first to relate it to another quan- 
tity associated with the same process X(t), but now in 
a semi-infinite geometry [0, oo]. Consider a process X(t) 
in [0, oo], starting at x and absorbed at the origin for 
the first time at tf— the first-passage time. Let m de- 
note the maximum of this process till tf (see Fig. [2]left). 
Then, it is clear that 1 — Q(x,L), the probability that 
the particle exits the box through the origin (and not 
through L), is precisely equal to the probability that the 
maximum to of the process in [0,oo] till tf stays below 
L, i.e., the cumulative distribution of to, Prob[m < L\x], 
in the semi-infinite geometry. The distribution of to is, 
in turn, related to the distribution of the first-passage 
time tf. Let q(x,T) = Prob[ij > T\x] denote the cumu- 
lative probability of tf, which is also the survival prob- 
ability of the particle starting at x in the semi-infinite 
geometry. One knows that for generic self-affine pro- 
cesses q(x,T) = q{x/T H ). For large T, q(x,T) ~ T~ e , 
where 8 is the persistence exponent of the process [T2] , 
This implies the scaling function q(y) ~ y e / H for small 
y [7J. Noting that to ~ t? for self-affine processes, it fol- 
lows that Q(x,L) = 1 — Prob[TO < L\x] = Prob[m > 
L\x] w Prob[i/ > L 1 '"^] = q[x/L]. This demon- 
strates the scaling behavior anticipated before, namely, 
Q(x,L) = Q(x/L), where Q(z) = q(z). Moreover, since 
q(y) ~ y e l R for small y, we get Q(z) ~ z^ for small 
z, with (f> = 8/H. For example, for Brownian motion 
H = 1/2 and 8 — 1/2, hence <f> = 1, in accordance with 
the exact result Q(z) = z. For the subclass of self-affine 
processes with stationary increments, the same exponent 
4> happens to describe the vanishing of the probability 
density close to an absorbing boundary [7J. 

Our general prediction Q(z) ~ z^ for small z can 
be verified explicitly for some self-affine processes where 
Q(z) can be computed exactly, as discussed later. More- 
over, we have numerically verified that this conjecture 
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FIG. 3: Behavior of Q(z) close to z = for fBm processes. 
For H = 2/3 (<j> = 1/2) the size of the box is L = 50, 200; 
for H = 3/4 {4> = 1/3) the size of the box is L = 100, 300. 
The continuum limit is reached when L — + oo. The expected 
slopes are reported as solid lines. Data have been shifted to 
make visualization easier. 



holds also for the fractional Brownian motion (fBm), i.e., 
a self-affine Gaussian process defined by the following au- 
tocorrelation function 



(t{ H + tl H 
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with < H < 1 [T3J [T3] . In [7J , we have proposed fBm 
as a natural candidate for describing the time evolution 
of the translocation coordinate. For this process, the 
persistence exponent is known, 8 = 1 — H [TJ], so that 
(j> = (1 — H)/H. An expedient algorithm for generating 
fBm paths is provided in |15j . The probability Q(z) can 
be numerically computed as follows. Given a realization 
of the process starting from the origin, we record its 
minimum and maximum values for increasing time; the 
process is halted when A" max — AT m i n > L. If the last up- 
dated quantity is X m ; n , the contribution to Q(x, L) is 
for x G (0, L — X max ) and 1 for x € (L — X max , L) . In the 
opposite case, the contribution is for x € (0, — X m - ln ) 
and 1 othewise. All simulations are performed by 
averaging over 10 6 samples. Fig. [3] shows the agreement 
between numerical simulations and predicted scaling of 
Q(z) for different values of the parameter H. 

Disordered potential. We next consider the stochastic 
motion of a single particle diffusing in a potential V(X), 
starting at the initial position x. The dynamics is gov- 
erned by the Langevin equation X(t) = f[X(t)] + r)(t), 
where X(0) = x and f(X) = -dV(X)/dX is the force 
and t](t) is a Gaussian white noise with (t](t)) = and 
(r](t)rj(t')) = S(t — t'). To compute Q(x,L), the first step 
is to write a differential equation satisfied by Q(x, L), tak- 
ing the initial position s as a variable and keeping the box 
size L fixed. Let us consider a small time interval [0, At] 
at the beginning of the process. In this time interval, the 
particle moves from its initial position x to a new position 
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FIG. 4: The difference between Eq. (pi and simulated Q(z). 
For fBm processes: H = 2/3 ((f) = l/2)with box size L = 200, 
and # = 3/4 (0 = 1/3), with box size L = 300. For fBm 
disordered potentials: Hy = 2/3 (<f> = 1/3), with box size 
L — 10 4 . For comparison, we display also the Sinai model 
H v = 1/2 {4> = 1/2), with box size L = 10 4 . 



x + Ax at time At, where Ax = f(x)At + r)(0)At, ?y(0) 
being the noise variable that kicks in at time 0. Since the 
process is Markovian, the subsequent evolution does not 
know about the interval [0, At], hence one gets 



Q(x, L) = (Q(x + f(x)At + T)(0)At, L)), 



(2) 



where () denotes the average over the initial noise ?7(0). 
Expanding the rhs of Eq. Q as a Taylor series in powers 
of At, using (r?(0)) = and (ry 2 (0)) = 1/At (delta cor- 
related noise), yields an ordinary differential equation, 
\Q"(x) + f(x)Q'(x) — 0. Solving with boundary condi- 
tions Q(0, L) = and Q(L, L) — 1 gives the exact result 



Q(x,L) 



J^e 2V ( x ')dx r 
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valid for arbitrary potential V(X). Note that for a 
potential- free particle, i.e., V(X) = in Eq. ([3]), we re- 
cover the Brownian result, Q(x, L) — x/L. 
Taking derivative with respect to x gives 
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2V(x) 



Peq ( X ,L) = -Q( x ,L) = ^ e2v{xi)dx . 



(4) 



which can be interpreted as the equilibrium probability 
density of the particle to be at x in presence of a poten- 
tial —V(X). When V(X) is a realization of a disordered 
potential, it is natural to introduce Q(x, L), the disorder- 
averaged hitting probability. An example where we can 
determine Q(x, L) explicitly is the classical Sinai model, 
i.e., when the potential V(X) is a trajectory of a Brow- 
nian motion in space, V(X) ~ A 1 / 2 \Ijg\. For this model 
the p oq (x,L) can be computed exactly [TTIHB] 



7T y/ X {L - X) ' 
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Thus, Q(x,L) = Q(z = x/L) again satisfies the generic 
scaling with a scaling function form with 



Q(z) = — an -in 

7T 
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Note that close to the origin Q(z) ~ with <f>— 1/2. On 
the other hand, it is well known that in the Sinai potential 
the particle evolves very slowly with time, X ~ hi 2 (t), 
showing a self-affine scaling in the variable T = logt, 
with H = 2. For this model, it is also known that the 
survival probability decays as 1/logt, i.e., T~ 6 , with 9 = 
l[19j[20]. Thus, 9/H = 1/2 = <p, in accordance with our 
general scaling prediction. 

We next consider a generic self-affine potential, 
V(X) - X Hv (with V(0) = 0), the Sinai potential 
being a special case with Hy = 1/2. We show that 
p eq (x,L) for such a potential is related to the proba- 
bility density of the location x m of the maximum of 
the potential V(X) over X e [0, L]. We rewrite Eq. 
Q as p cq (x,L) = [J L e^^'l-v^dx'}- 1 , rescale vari- 
ables x' — > x' L and x — > xL and use the self-affine 
property V(xL) — L Hv V(x) to obtain p cq (x,L) = 
[ji e 2L H v(v(x')-v(x)) dx i}-i_ large L, using a steepest 
decent method, we immediately sec that, for each realiza- 
tion of the disorder potential V(X), p cq (x, L) ~ S(x—x m ) 
where x m denotes the position where V(X) is maxi- 
mum. This observation has two immediate consequences: 

(i) By integrating over x, we get, for each realization, 
Q(x,L) ~ 9{x — x m ). This means that, for any given 
realization, if the starting position x is to the left (right) 
of the location x m of the maximum, Q(x, L) ~ (respec- 
tively Q(x, L) — > 1) indicating that the particle exits the 
box through (through L) as depicted in Fig. fright). 

(ii) By taking average over the disorder, we get 



p eq (x,L) ~p m (x,L) 



(7) 



where p m (x,L) is the probability density that the max- 
imum of the potential V(X) over [0, L] is located at x. 
For example, for the Sinai case, one knows from Levy's 
arcsine law [5] that p m (x,L) = 1 fityj x(L — x). Thus, 
in this case the relation ^ is verified by the exact re- 
sult (|5]). However, the relation ([T]) is more general and 
holds for arbitrary self-affine potential. We remark here 
that physically the relation ^ reflects the fact that in 
a self-affine potential —V(X) in which the particle is at 
equilirium, the limit L —>■ oo limit is equivalent to the 
zero temperature T — > limit forcing the particle to the 
minimuum of the potential —V(x) or equivalently to the 
maximum x m of V(x). 

A useful consequence of ^ is that it allows us to relate 
the persistence or the survival probability of a particle 
moving in a disordered self-affine potential to the sta- 
tistical properties of the potential V(X) itself. The dis- 
ordered potential V(X) ~ X Hv (we assume V(0) = 0) 
can itself be regarded as a stochastic process with the 
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space coordinate X playing the role of 'time'. So, the 
probability that V(X) stays below (or above) the level 
X = up to a distance L decays, for large L, as L~ 6v 
where 9y is the spatial persistence exponent [2T] of V(X). 
For example, for the Sinai potential (brownian motion in 
space), Qy — 1/2. The pair of exponents (Hy,9y) asso- 
ciated with the potential can now be related to the corre- 
sponding exponents associated with the temporal motion 
of the particle in this potential. By Arrhenius' law for 
the activated dynamics, the time required for particle 
diffusing in V(X) to overcome an energy barrier scales 
as t ~ e v(x \ Using V(X) ~ X Hv , one deduces that 
X ~ T x l Hv where T = log(t). Thus the particle mo- 
tion X(T) ~ T H is a self-affine process as a function of 
T = log(i), with a Hurst exponent H = 1/Hy. Next, 
we note that p m (x,L), the probability that the maxi- 
mum of V(X) occurs at x, coincides, when x — > 0, with 
the probability that V(X) stays below up to a dis- 
tance L, hence p m (x — > 0,L) cx L~ 0V . On the other 
hand, based on our general argument, we expect that 
Q(x,L) ~ Q/L)^ wh en x -> 0, where = 0/if. this 
means that p eq (x,L) cx x^ 1 !^. Note that here 6* is 
the persistence exponent associated with the temporal 
motion of the particle, i.e., the survival probability of 
the particle up to time T = log(i) decays as ~ T~ e . 
Matching powers of L from both sides of ^ provides the 
desired relation between temporal and spatial exponents 
9 = OyH = 8y /Hy. For instance, in the Sinai model, us- 
ing 9y = 1/2, Hy = 1/2 we get and 9 = 1, in agreement 
with the exact result [19, 20 . If one considers a potential 
V"(X) = £(X) where £(X) is a white noise in space, it 
is self-affine with Hy = 3/2. The exponent 9y = 1/4 is 
known exactly |26j . Thus we predict that for this poen- 
tial, the survival probabaility up to time t will decay as 
- {\ogt)- e with = 1/6. 

Super-universality of Q(z). For some non-Brownian 
stochastic self-affine processes, the full function Q(z) is 
known. For instance, Levy Flights are Markovian su- 
perdiffusive processes whose increments obey a Levy sta- 
ble (symmetric) law of index < /i < 2. The Hurst 
exponent is H = l//i. By virtue of the Sparre Andersen 
theorem [55] , the persistence exponent is 9 — 1/2, inde- 
pendent of /i. Hence, <f> = 9/H = fi/2 (see also [2"51 l2"4"] j. 
The full function Q(z) for Levy Flights has been com- 
puted [25 and can be recast in an elegant form 

Q(z) = 4(0, <t>) = P|| j" Ml - u)f- x du, (8) 

i.e., a regularized incomplete Beta function containing 
a single parameter 4> = /x/2. Clearly, Q(z) ~ z^ as 
z — > in agreement with our prediction. The formu- 
lae for Brownian motion (with (j) — 1) Q(z) = z and 
for the Sinai model {cf> = 1/2) in ^ can also be ex- 
pressed as Moreover, the distribution of the maxima 
for a symmetric Levy Flight process is given by (TsJ) , by 
virtue of the Sparre Andersen theorem |22j . Hence, we 



expect the hitting probability ^ to apply also to parti- 
cles diffusing in a Levy Flight disordered potential, with 
4> = 1/2. Finally, Q(z) is known also for the Random 
Acceleration model, a non-Markovian process that is de- 
fined by d 2 X/dt 2 — r/(t), with rj(t) as before. The motion 
starts at X(0) = x, with initial velocity v(0) — 0, and 
is superdiffusive, with X ~ i 3 / 2 , i.e., H — 3/2. Its first- 
passage properties have been widely studied [26 . The 
persistence exponent is 9 = 1/4, so that <f) = 9/H = 1/6. 
Bicout and Burkhardt [57] also computed the full exit 
probability Q(x,L). One can again recast this formula 
in the same super-universal form ^ with = 1/6. 

Based on these special cases, it may be tempting 
to conjecture that the full function Q(z) for arbitrary 
anomalous doffusion processes has the super-universal 
form ([8]), the only information about the process enters 
this formula through a single exponent <p. However, this 
turns out not to be the case, and we are able to show no- 
table counterexamples. In Fig. [4] we compute the hitting 
probability for fBm self-affine processes and for particles 
diffusing in fBm disordered potentials, and display the 
numerical difference with respect to formula ([8]), with 
the appropriate exponent <f>. We find that in neither case 
Q(z) can be described by the super-universal form ([8]). 
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